Current-driven vortex oscillations in metallic nanocontacts 



Q. Mistral, 1 M. van Kampen, 2 G. Hrkac, 3 Joo-Von Kim, 1 ^ T. 

Devolder, 1 P. Crozat, 1 C. Chappert, 1 L. Lagae, 2 and T. Schrefl 3 

1 Institut d'Electronique Fondamentale, UMR CNRS 8622, Universite Paris-Sud, 91405 Orsay, France 
2 IMEC, Kapeldreef 75, B-3001 Leuven, Belgium 
J Department of Engineering Materials, University of Sheffield, Sheffield SI 4DU, United Kingdom 

We present experimental evidence of sub-GHz spin-transfer oscillations in metallic nano-contacts 
that are due to the translational motion of a magnetic vortex. The vortex is shown to execute 
large-amplitude orbital motion outside the contact region. Good agreement with analytical theory 
and micromagnetics simulations is found. 
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Lateral confinement in magnetic nanostructures leads 
to the appearance of complex magnetic states. One 
pertinent example concerns the magnetic vortex, which 
is well-studied and has gained renewed interest due to 
improved lithography and fabrication techniques. The 
dynamics of such topological objects is, in part, gov- 
erned by the confining potential. For instance, mag- 
netic vortices possess a translational mode by which 
the vortex core executes a spiralling motion about its 
equilibrium position [TJ [2]. Recent experimental stud- 
ies have demonstrated resonance phenomena with such 
cores [31 HI 03 O [7], and proposals have been made for 
magnetic resonators based on these objects. 

More recently, it has been shown experimentally that 
such vortices can be brought into a self-oscillatory 
state [5]. This has been achieved with the spin-transfer 
effect in which a spin-polarized current transfers spin- 
angular momentum to the magnetic structure, giving 
rise to an additional torque exerted on the magnetiza- 
tion j9[ [10]. Under certain conditions the spin-transfer 
compensates the natural dissipation processes, leading 
to self-sustained oscillations. In contrast to a resonator 
in which oscillations are driven by a periodic external 
force, the self-oscillator here is driven by a constant 
electric current only. While current-driven oscillations 
related to spin-waves in magnetic multilayers are well- 
studied [TTJ [TSJ [T3] , Ref. [5] provides the first conclusive 
evidence of this phenomenon for a magnetic vortex. 

In the context of spin-transfer, confinement can also 
arise from the distribution of applied current itself. One 
pertinent example is the metallic point-contact in which 
the current is applied through a small metallic cross- 
section in contact with a continuous magnetic film |13j . 
In the absence of any applied currents, a vortex can exist 
in a continuous film but, aside from sample edges and 
point defects, there is no magnetic potential in which 
vortex oscillations can take place. In the presence of 
an applied current, however, the Oersted fields produce 
an attractive potential centered at the contact center. 
Earlier experimental studies have provided hints at the 
existence of vortex modes, but no explanations for the 



observed low-frequency excitations were given |14j . In 
this Letter, we show that the cylindrical nature of the 
Oersted fields, generated by the applied current, is favor- 
able to the creation of a magnetic vortex and provides 
a potential in which the vortex oscillates. Moreover, we 
demonstrate that the sub-GHz frequencies of these vortex 
self-oscillations result from a vortex orbit that is actually 
outside the contact region. 

The magnetic structure we have studied in experiment 
is a metallic point contact deposited on a metallic spin- 
valve stack. The multilayer is composed of Ta(3.5 nm) / 
Cu(40) / Ta(3.5) / Ni 80 Fe 20 (3) / IrMn(7) / Co 90 Fe 10 (3.5) 
/ Cu(3) / Ni 80 Fe 20 (4) / Pt(3). The NiFe (4 nm) layer 
is the magnetic free layer and the CoFe layer, which is 
exchange biased by the IrMn antiferromagnet, serves as a 
reference layer for the giant magnetoresistance and spin- 
transfer. The layers are sputter-grown in an ultrahigh 
vacuum system with a base pressure of 3 x 10 -8 Torr. 
The metallic point contacts are made to the laterally ex- 
tended (> 5 /im) stack which is incorporated in the cen- 
tral conductor of a coplanar waveguide (CPW) structure. 
After patterning by conventional lift-off lithography, the 
stack is passivated by a 50 nm Si02 layer deposited using 
rf sputtering. The point contacts are defined in a PMMA 
layer by electron-beam lithography and etched into the 
Si02 by a short dip in a buffered HF solution. Next, 
the CPW structure is deposited and patterned by lift-off 
lithography. Finally, the devices are annealed at 250 °C 
for 10 minutes to improve the exchange bias coupling 
between the IrMn and the CogoFeio layers. 

The properties of the point contact have been charac- 
terized in detail by scanning electron microscopy (SEM), 
atomic force microscopy (AFM), ellipsometry and opti- 
cal microscopy measurements. The target diameter of 
the point contacts are nominally 100 nm, however, due 
to the isotropic wet etch the diameter of the contacts is 
increased with respect to the 100 nm holes in the PMMA 
layer. Indeed, our measurements show that there is a sig- 
nificant variation in the diameter as a function of depth. 
In Fig. [pi and b, we present an SEM and AFM measure- 
ment, respectively, of a typical point contact. While a 



diameter of 232 nm is extracted from SEM, AFM mea- 
surements show that the diameter can vary from 130 to 
270 nm over a thickness of 40 nm. Based on the depth 
profile and our multilayer configuration, we estimate the 
contact diameter at the free layer to be approximately 
200 nm. 

For the magnetic properties, a uniaxial anisotropy field 
of /Zo-Hfc — 1.4 mT for the free layer is determined from 
hard-axis hysteresis loop measurements. From the easy 
axis loops, we determine the biased layer switching to be 
around —60 mT. A small bias of 6 mT in the free layer 
switching is also observed and is attributed to Neel cou- 
pling between the magnetic layers. The saturation mag- 
netization of the free layer is found to be HqM s = 1.1 T 
from alternating gradient field magnetometer measure- 
ments. For the electrical measurements, the applied cur- 
rent perpendicular through the spin- valve stack is gen- 
erated by applying a DC voltage to a 200 resistor in 
series with our sample. A 20 mfi magnetoresistance is 
a typical figure for our spin-valve stacks. Two bias tees 
allow DC and RF routing. One of the RF routes is am- 
plified (+43 dB over the 100 MHz - 26 GHz band) and 
connected to a spectrum analyzer. The other RF route 
is terminated by a 50 Ohm load while the DC route is 
shorted. The final power spectra data is obtained after 
subtracting a reference curve taken with a zero current 
bias. The rf measurements are made for magnetic fields 
applied perpendicular to the film plane, which improves 
the magnetoresistance signal because the projection of 
the free layer oscillations onto the fixed layer is larger. 
Positive currents describe electrons flowing from the fixed 
to the free layer. 

An example of low-frequency power spectra is pre- 
sented in Fig. [TJ;. The peaks shown represent the fun- 
damental mode for several applied currents under a per- 
pendicular field of noH — 210 mT. Harmonics to fourth 
order are also observed, but these are several orders of 
magnitude smaller than the fundamental, which indicates 
magnetization precession with small ellipticity. The oscil- 
lations have quality factors up to 250, which corresponds 
to the narrowest linewidth of 1.2 MHz observed at a fre- 
quency of 299 MHz. These low precession frequencies 
are suggestive of vortex motion rather than uniform pre- 
cession [6j. Indeed, the Smit-Beljers equation fTS] for 
uniform magnetization precession predicts a precession 
frequency of 1.08 GHz with our experimental parame- 
ters, which is clearly four to five times higher than what 
we observe. In contrast to earlier experiments in which 
nominally smaller diameter contacts were studied |13j . 
the Oersted fields in our samples are significantly larger 
(at equivalent current densities) at the contact edges and 
therefore favor vortex formation. For example, a 30 mA 
current generates an Oersted field of 60 mT at the contact 
edge, which is a significant fraction of the perpendicular 
applied fields we considered. 

Two other key experimental observations support the 



2 




Frequency (GHz) 

FIG. 1: (Color online) (a) SEM and (b) AFM image of 
point contact, (c) Low frequency power spectra measured 
at hqH = 0.21 T for several applied currents. The field and 
film geometry is shown in the inset of (c). 



hypothesis of vortex oscillations. First, the measured os- 
cillation frequencies increase as a function of applied cur- 
rent for all applied fields considered (hqH < 350 mT). An 
example of this behavior is presented in Fig. [2] in which 
the full power spectra is shown in a density plot. This fre- 
quency blucshift is a feature of magnetization precession 
out of the film plane, but is not expected for the fields 
used in experiment. Indeed, the spin-torque oscillator 
theory of Ref. 1161 for uniform magnetization precession, 
predicts a frequency redshift for these applied perpen- 
dicular fields. Second, the precession mode is sensitive 
to small in-plane applied fields. The oscillations are ob- 
served to be stable for applied in-plane fields up to 15-20 
mT in magnitude, and for larger fields the oscillations 
are completely suppressed. Furthermore, the narrowest 
peaks are obtained for in-plane fields that compensate 
the 6 mT bias field due to the Neel coupling, which is 
again consistent with the vortex picture. 

To understand the exact nature of the low frequency 
oscillations, we performed full micromagnetic simulations 
of a trilayer CoFe(3.5 nm) /Cu(3) /NiFe(4) stack. The sat- 
uration magnetization of the free layer is taken to be the 
same as the experimental value, /j,qM s = 1.1 T, and a 
GMR ratio of 1% is used. We account for the inhomoge- 
neous current distribution flowing through the free layer 
by computing the local current density from the local an- 
gle between the free and fixed layer magnetizations. The 
Oersted field is computed with the Biot-Savart law from 
this current distribution and an asymmetric Slonczewski 
term for the spin-transfer is used [17] • The simulation 
area is a circular film 1000 nm in diameter, with the 200 
nm diameter point contact at the center. The system 
is discretized with a finite element method with a linear 
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FIG. 2: (Color online) (a) Color map of experimental power 
spectral density (PSD) as a function of current for /j,oH = 
350 mT. Solid squares represent results of micromagnetics 
simulations, (b) Top view of vortex magnetization profile 
obtained from simulation, (c) In-plane component of magne- 
tization represented in a gray level (black for m y — +1, white 
for m y — — 1) for the entire simulation area. Also shown are 
the contact and the vortex orbit. 



basis function and a discretization size of 4 nm, which is 
below the exchange length of 4.5 nm. An exchange bias 
field of 162 mT is applied to the reference CoFe layer to 
simulate pinning in the experimental stacks. 

The simulations are performed as follows. First, we 
calculate the remanent state with an external field ap- 
plied perpendicular to the film plane and in the absence of 
currents. From an initial uniform state in the film plane, 
we allow the magnetization to relax toward the applied 
field direction through time-integration of the equations 
of motion. During this process, a vortex first forms under 
the point contact (which is situated in the center of the 
simulation grid) and is followed by the creation and anni- 
hilation of multiple vortices. After this transient phase, 
a single vortex state emerges as the system relaxes to its 
remanent state. Next, a current is applied through the 
point contact with currents ranging from 26 mA to 36 
mA, with the remanent state serving as the initial state 
for this step of the calculation. We observe that the ad- 
ditional spin torque drives the vortex out of the contact 
area and towards a stable orbit around the contact. 

These simulations show that the oscillations observed 
in experiment are related to a large-amplitude transla- 
tional motion of a magnetic vortex. In contrast to the 
vortex dynamics observed in nanopillars in which the vor- 
tex core precesses within the confining part of the Oer- 
sted field [5] , the dynamics here correspond to an orbital 



motion outside the contact region. This behavior can be 
likened to planetary orbital motion under the influence 
of a gravitational field; the spin-transfer torque leads to 
a centripetal motion of the vortex core, which is coun- 
terbalanced by the attractive potential provided by the 
Oersted field. Good quantitative agreement between the 
simulation and experimental frequencies is achieved, as 
shown in Fig. [2j In addition, the stability of the vortex 
has also been investigated using simulations. For small 
in-plane fields (less than 10 mT for this system), a small 
blueshift in the frequency is observed. This leads to an 
asymmetry in the field profile relative to the vortex struc- 
ture and results in an elliptical motion of the vortex orbit. 
We have also verified in simulation that larger in-plane 
fields lead to a breakdown of the vortex structure. 

A qualitative understanding of how the oscillation fre- 
quencies change with applied field and current can be 
achieved in terms of a rigid vortex model. Let 0(r) and 
$(r) denote the magnetization orientation in spherical 
coordinates, with the spatial variations in r given in po- 
lar coordinates (p, if) to exploit the symmetry of the cir- 
cular point contact. We take the magnetization to be 
uniform across the film thickness. The vortex profile is 
approximated by a polar angle that only depends on the 
radial spatial coordinate, O = O(p), and an azimuthal 
angle that depends only on the azimuthal spatial coordi- 
nate <& = ip ± 7r/2, where the sign in the last expression 
denotes the vortex polarization. 

Based on a given magnetization profile, we can derive 
an equation of motion for the vortex core by following 
Thicle's method [18], 
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where Xq = (Xq,Yq) represents the vortex core posi- 
tion in the film plane, 7 is the gyromagnetic ratio, a is 
the Gilbert damping constant, a is a measure of spin- 
transfer efficiency, G = J dVsin©(V<i> x V6) is the 
gyrovector, D = J dV (V0 VO + sin 2 OV$ (g) V$) 
is the damping dyadic, and W is the total vortex en- 
ergy. The new terms due to spin-transfer are Pm = 
P\\ Ipc dV ( V ® sin* + |V<& sin 20 cos $) , which is the 
component of spin-transfer parallel to the film plane, and 
P± = pj_ J pc dV sin 2 O V<&, which is the component of 
spin-transfer perpendicular to the film plane. The unit 
vector (p\\,p±) represents the orientation of the incom- 
ing spin-polarization and the subscript PC denotes an 
integral over the region underneath the point contact. 

In the absence of pinning centers and assuming the 
sample edges are sufficiently far from the vortex core, 
the only position-dependent component of the vortex 
energy arises from the Oersted field generated from 
the current flow through the point contact. Assum- 
ing the current flows through a perfectly cylindrical 
cross-section in the free layer, we obtain for W in 
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polar coordinates, W = — /irjM s /fj J dVf(f) sin O x 
(sin tp cos $ — cos (p sin <&) , where f(f) describes the ra- 
dial variation of the Oersted field amplitude, Hi, from 
the contact center. While it is difficult to evaluate this 
energy exactly, we find numerically that the functional 
form W(p) ~ nip gives a reasonably good approxima- 
tion for realistic O profiles. 

The vortex oscillation frequency is found by solving 
Eq. [T] in the rigid vortex approximation. Assuming the 
vortex core is sufficiently far from the contact region, we 
can suppose that the gradient in the polar angle van- 
ishes, VO = 0, and the polar angle is given by the uni- 
form tilted angle due to the perpendicular applied field, 
cos 0o = H/M s , within the contact region. Evaluating 
the integrals in within this approximation and using 
the representation i?o ex p(*V-') = Ao + iYo, we find the 
simplified equations of motion, 
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where G and D are constants corresponding to the gy- 
rovector and damping dyadic, respectively. We are inter- 
ested in the steady state precession for which Rq — and 
u> = ip is a constant. Using these constraints, we find a 
vortex precession frequency of 

2 aDn 2 \ J 

—O 9 • ( 4 ) 

ja z crp± J sin O 

The dependence on the applied field strength can be in- 
troduced by substituting 9 = cos -1 (H/M s ) and approx- 
imating the perpendicular spin-polarization component 
from the z component of the pinned-layer magnetization, 
p± ~ H/M2, where M2 is the saturation magnetization 
of the pinned layer. Note that only the perpendicular 
component of the spin-polarization leads to vortex pre- 
cession; the in-plane component only leads to a displace- 
ment of the vortex core. In our system, a non-vanishing 
contribution from p± arises from the finite tilt angle of 
the fixed layer magnetization by virtue of the perpendic- 
ular applied field. 

A comparison between theory and experiment is shown 
in Fig. [3] As the current is varied, we observe a quasi- 
linear variation of the oscillation frequency. This behav- 
ior is consistent with theory and is in essence due to the 
linear dependence of the Oersted field energy on current. 
To produce the theoretical curves in Fig. [3^,, we first fit 
the H = 105 mT series using Eq.|4]and then use the theo- 
retical H dependence to compute the slopes for the other 
applied fields, which are in good quantitative agreement 
with the experimental data. The theory also accounts 
for the frequency variation with applied field, as shown 
in Fig. [3J5. As before, we fit Q to one data set (I = 
36 mA) and used the predicted current variation for the 
other sets. 
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FIG. 3: (Color online) Vortex oscillation frequency as a func- 
tion of (a) applied current and (b) applied field. Dots repre- 
sent experimental data and solid lines are fits based on Eq. [4] 



In summary, we have shown with experiment, sim- 
ulation, and analytical theory that current-driven sub- 
GHz oscillations in point-contacts are a result of large- 
amplitude vortex dynamics. The vortex attains a stable 
circular orbit outside of the contact region, whose fre- 
quency is governed by the applied perpendicular field and 
current. 
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